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Abstract 

We develop the formalism for the one-loop no-boundary state in a cos- 
mological model with fermions. We use it to calculate the reduced density 
matrix for an inflaton field by tracing out the fermionic degrees of freedom, 
yielding both the fermionic effective action and the standard decoherence 
factor. We show that dimensional regularisation of ultraviolet divergences 
would lead to an inconsistent density matrix. Suppression of these diver- 
gences to zero is instead performed through a nonlocal Bogoliubov trans- 
formation of the fermionic variables, which leads to a consistent density 
matrix. The resulting degree of decoherence is less than in the case of 
bosonic fields. 
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1. Introduction 



Quantum cosmology - the application of quantum theory to the Universe as 
a whole - is undergoing now a period of successful development connected with 
the fact that it can help to explain the origin of inflation in the early Universe. 
One of the basic problems in quantum cosmology is the choice of an appropriate 
boundary condition. Two proposals have become especially popular: The no- 
boundary state m, and the tunneling state || . Analysis of these two proposals 
in the one-loop approximation allows one to formulate the conditions for the 
normalisability of the wave function of the Universe ||J] and to calculate the 
probability peak for the energy scale of inflation in different models || |6| . Such a 
probability peak has been interpreted as providing a criterium to select amongst 
the members of an "ensemble" of classical universes. 

However, quantum theory does not yield a classical ensemble. A complete 
analysis would thus have to include a quantitative discussion of the quantum-to- 
classical transition. 

It is now generally accepted that the classical properties for a subsystem 
arise from the irreversible interaction of this quantum system with its natural 
environment. Starting with the pioneering work of Zeh in the seventies [[T|, this 
concept of decoherence has been developed extensively, see [|J and [§] for reviews. 
Quite recently, this continuous loss of coherence was observed in quantum-optical 
experiments [|10] . 



Quantum-to-classical transition through decoherence also applies to quantum 
cosmology where it was suggested to consider background degrees of freedom 
(such as the scale factor and the homogeneous mode of the inflaton scalar field) as 
the relevant system and the perturbations (such a gravitational waves or density 
fluctuations) as an environment |Tl| . A quantitative discussion of this idea was 
first done in El , where it was demonstrated how the scale factor and the inflaton 
field can acquire classical properties. This idea was further pursued in many 
papers M. The quantum-to-classical transition through decoherence plays also a 
crucial role for the primordial fluctuations itself that eventually serve as classical 
seeds for galaxies and clusters of galaxies . 

However, there is one big problem arising in the application of this idea to 
cosmology. Since there are infinitely many environmental modes, there arise ul- 
traviolet divergences that have to be removed. It was already noted in |14[] that 
dimensional regularisation is only applicable for the phase part of decoherence 
factor, not its absolute value. In our preceding paper |15j this was confirmed 



by detailed calculations. It was also shown in [p~5|| , that for the case of bosonic 
fields, using the so-called conformal parametrisation one can get rid of ultraviolet 
divergences and obtain finite results for the absolute value of the decoherence fac- 
tor. An analogous idea was investigated in QED and in a particular cosmological 



model in |L6| . A similar proposal was also discussed in [T] 

The purpose of the present paper is to extend the analysis of [I5| to fermionic 
degrees of freedom. From a physical point of view, this is important because 
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fermions constitute the main part of physical fields in our Universe. On the 
other hand, the Grassmannian character of fermion variables implies an essential 
modification of the formalism which is of interest by itself. 

First of all, we want to emphasise that the formalism of the Hartle-Hawking 
wave function of the Universe via a path integral in one-loop approximation was 
previously developed only for bosonic fields [|], |, Here, we extend this for- 
malism to fermionic fields. In this case the structure of the norm of partial wave 
functions, describing perturbations of different harmonics of fermionic field - an 
important ingredient of the one-loop formalism - is more complicated than in the 
case of bosonic fields where it can be reduced to calculating the Wronskian of 
the corresponding basis functions. Using this formalism, we can then calculate 
the decoherence factor for the reduced density matrix of scale factor and inflaton 
field. However, like in the bosonic case, ultraviolet divergences appear, as was 



already noted in an early paper on decoherence due to fermions [[18 



we 



By handling the divergent series within dimensional regularisation [ F9 , 2"H 
obtain a renormalised expression for the absolute value of the decoherent factor. 
Unfortunately, however, it shows a wrong behaviour violating the consistency 
conditions for a density matrix as a bounded operator, let alone its decoherence 
properties. A similar phenomenon was observed also in the case of bosonic fields 



15[| . However, a simple reparametrisation of field variables by overall multipli- 



cation with an appropriate power of the scale factor allowed us there to obtain 
finite expressions. This does not work in the case of fermionic fields, where 
such a reparametrisation is incompatible with the structure of the measure on 
the corresponding Hilbert space of states. However, there exists another pos- 
sibility: Performing a Bogoliubov transformation of field variables one can get 
again finite expressions for the reduced density matrix. It is shown that such 
Bogoliubov transformations are in a certain sense analogous to the well-known 
Foldy-Wouthuysen transformation of spinor variables [2I|. Thus, as in the case 



of bosonic fields, the definition of the environment is crucial for the consistency 
and decoherence properties of the reduced density matrix, but in the fermionic 
case this definition is much more involved. 

Our paper is organised as follows. In section 2 we briefly review the spinor 
formalism for quantum cosmology. This section includes the relation between the 
two- and four- component spinor formalism and presents the expansion of Weyl 
spinors in three-dimensional spinor harmonics with a special emphasis on the 
dimensionalities of the 0(d) spinor representation on ^-dimensional spheres (nec- 
essary for dimensional regularisation). Section 3 is devoted to the quantisation 
of fermion fields in the Hamiltonian formalism and to the construction of the 
cosmological wave function. In section 4 we discuss the connection between path 
integral, Euclidean effective action and density matrix of fermions. Section 5 con- 
tains the presentation of contributions of massless and massive fermions into the 
amplitude of the decoherence factor in the reduced density matrix. In section 6 
we demonstrate that the application of dimensional regularisation violates basic 
properties of the density matrix. Section 7 comprises the discussion of Bogoliubov 
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transformation of fermionic variables, allowing to get rid of ultraviolet divergences 
in the decoherence factor. We also exhibit the connection with the traditional 
Foldy-Wouthuysen transformations. Section 8 contains a brief summary of the 
obtained results. 



2. Spinor formalism for quantum cosmology 

In the main part of this paper we shall use a two-component spinor formalism 



as it is usually done in papers with applications to quantum cosmology [22, 23. 24 



However, at the beginning we would like to present some formulae connecting this 













= Mf^ B 







formalism with that of four- dimensional Dirac spinors (see, for example Ref. |25fl). 

We shall use two-component Weyl spinors in Van der Waerden notation. 
There are two types of spinors with primed and unprimed indices which are 
transformed under the action of matrix M £ SL(2, C) according to the following 
rules: 



(2.1) 

Those with dotted indices transform under the (o, |) representation of the Lorentz 

group, while those with undotted indices transform under the f§,0j conjugate 
representation. For the Minkowski metric we use the signature 

r] mn = diag(-l, 1, 1, 1). (2.2) 

The Pauli matrices a^s have the following form: 

Indices can be lowered and raised by tensors s A b an d e AB with components 
£il = ^22 = £ n = £ 22 = 0, 

e 21 = e 12 = l, e 12 = e 21 = -l. (2.4) 

(The components of e with dotted indices are the same). It is convenient also to 
define 

-mAA = £ AB £ AB ^m. (25) 

The anticommutation relations for the Pauli matrices are: 
(a a +a a ) A = -2r] d A , 

[a a + cr a ) B = —2rj d B . (2.b) 
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For four- dimensional 7 - matrices and Dirac spinors we shall use the Weyl 
basis. In this basis, 7 - matrices are given 



7™ = ^„ j , (2.7) 
while the Dirac spinor is related to the two-component Weyl spinors by 

* = ( X A) ■ (2-8) 




It is worth noting that the canonical basis for the Dirac 7 - matrices, 

■ 0) ■ M 

is related to the Weyl basis (|2.7| ) by the following similarity transformation: 

T w = xv c x~\ x = ^= I 1 ']) . (2.10) 



Note that the Weyl basis is more consistent from the group-theoretical point of 
view because in this basis the Dirac spinor is represented as a direct sum of two 
spinors belonging to irreducible representations of the Lorentz group (or to be 
precise, its covering group SL(2,C)). This basis is also more convenient for the 
description of the massless limit for ultrarelativistic fermions. At the same time, 
the canonical basis is convenient for studying the non-relativistic limit and for the 
transition from the Dirac equation to the Pauli equation. We shall come back to 
this question when considering the relation between Bogoliubov transformations 
and the Foldy-Wouthuysen representation for Dirac spinors. 

Now, starting with the traditional action for Dirac spinors in four-dimensional 



spacetime [26, 21 1 and using the above relations (|2.7p,([2.8|), one can get the fol- 



lowing action | 22| : 



S = d^xei^e^D^+x^D^+h.c.-m J d A xe( X A<j> A +0 A XA), (2.11) 

where is a spinorial covariant derivative, : is the Pauli matrix contracted 
with the tetrad e(i, and e is the determinant of the tetrad. 



Varying the action (|2.11j) , one can get the Dirac equations 





= imXAi 






= im&Ai 






= -imxA, 




<A D »X A 


= —im(f)A- 


(2.12) 
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In what follows we shall consider fermionic perturbations on the background 
of a spatially closed Friedmann model with S 3 spatial sections. We shall work 
within the inflationary model with the spacetime metric || 

ds 2 = -dt 2 + a 2 (t)dQ 2 3 , (2.13) 

(dQl is the metric on the unit three-sphere) with a cosmological scale factor 
evolving according to the law 

a(t) = — cosh if*, (2.14) 
H 

where H = H(<p) is an effective Hubble constant possibly generated by an inflaton 
scalar field (p. 

Then one should expand the fermionic field using the full set of spinor har- 
monics on the three-sphere |22| : 



a" 3 / 2 



2tt 

^" np q 

a 3 / 2 



E E <Vn P 0)PA fa) + r np (t)al q (x)}, 



E E ocl q [rhn P {t)p n Hx) + r np (t)a^(x)}, 



A 27T 

" n np q 

a- 3 / 2 

XA = ^— EE^MOPIV) +t np (t)a n A %x)}, 

* np q 

3/2 

* np q 



Here, the weight factor a 3//2 is especially chosen in order to write the action ( 2.11 ) 



in a special parametrisation reflecting the conformal properties of the spinor field. 
It is conformally invariant in the massless case and, therefore, decouples in this 
parametrisation from the gravitational background. Summation over n runs from 
1 to 00, while the indices p and q run from 1 to n(n + 1) - the dimensionality of 
the irreducible representation for Weyl spinors. The time-dependent coefficients 
mnp,r n p,t n p, s np and their complex conjugates are taken to be odd elements of 
a Grassmann algebra. The constant coefficients are included for convenience in 
order to avoid couplings between different values of p in the expansion of the 
action. They may be regarded as block-diagonal matrices of dimension n(n + 1) 



whose structure can be found in ||22|| . The harmonics p,a,p,a are normalised 
eigenfunctions of the Dirac operator on the three-sphere with eigenvalues ±(n + 
1/2), n > 1. Here we are using a notation that differs from [^] by a unit shift 



of the quantum number n. Inserting ( |2.15| ) into the action ( |2.11| ) and using the 



properties of spinor harmonics, one finds that the action reads 



S = Y.jdt 



% , _ ± _ . 

2 v^wp^'Wp W^npWlnp SnpSnp ^np^np 



np 

~~\~t n pt n p ~\- t n pt n p -\~ T n pT n p ~\- T n pT n p} 
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~\~ [Jl ~\~ ) (jTlnp^Tlnp ^np^np ^np^np ^np^np) 



1, 1, 

—in + - 

Tn\T n pt n p -\- t n pT n p -\- S n pTfl n p -\- rfl n pS n p) 

~)~7) ^ ] (j^np^^np SnpSnp ^nptnp ^np^np) g, 
np 

^ ] (j^np^^np SnpSnp tnptnp T'npT'np) • (2-16) 
* np 

The last two lines in this equation represent the boundary terms at the initial (Si) 
and final (Sp) hypersurfaces of the spacetime domain. They have to be necessarily 
included in order to guarantee the consistency of the variational problem for this 
action with fixed unbarred variables at Sp and fixed barred variables at Si - the 
choice of boundary conditions characteristic for the definition of the fermionic 
wave function in the holomorphic representation [F22J. Thus the full action has 
the form of a sum of fermionic actions 

S n = J dt (^-{xx + xx + yy + yy) + e~ a (n + -){xx + yy) — m(yx + xy)j 

+-{xx + yy) Sp + -(xx + yy)s r (2.17) 

with the Grassmanian variables x and y denoting m nv and s np or t np and r np , 
respectively. 

From the action ( 2.17[ ) the following field equations may be derived: 

ix + vx — my = 0, (2-18) 
ix — vx + my = 0, (2-19) 
%y + vy + mx = 0, (2.20) 
iy — vy — mx = 0, (2.21) 

where 

n + 1/2 

v = 1 . (2.22) 
a 

In the rest of this section we give the formula for the dimensionality of 0(4) 
irreducible representation for spinors on spheres of arbitrary even dimensionality, 
which we shall need in the dimensional-regularisation calculations below. For 
this we use the technique which for scalar, vector and tensor harmonics on S 3 
has been constructed by Lifshitz and Khalatnikov [27 and for spinor harmonics 
on S 3 is described in |22j. 

Let us consider flat Euclidean four-space with the metric 

ds 2 = 6^dx"dx u . (2.23) 
The metric may also be written in spherical coordinates r, x, &, </> : 

ds 2 = dr 2 + r 2 dn 2 3 , (2.24) 
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where 



dnj = d X 2 + sin 2 X (d6 2 + sin 2 6d<p 2 ). (2.25) 

Harmonics on S 3 are then obtained by considering the following homogeneous 
polynomials in Cartesian coordinates, 

r n p A (x, 0, 0) = T AAi ... An _ iAi ... An _ iX A ^ ■ ■ ■ x ^-^-i , „ = !,..., (2.26) 

where 

X AA = e ^V (2.27) 

and T AA A _ i _ A is a constant spinor of rank (2n — 1), totally symmetric 
in all its indices and hence traceless (in the metric Eab)- It is obvious that the 
number of independent harmonics is defined by the number of different choices of 
dotted and and undotted indices which can acquire only two values: or 1. For 
n — 1 dotted indices there are n opportunities while for n undotted indices there 
are (n+1) different choices. Thus, the dimensionality of irreducible representation 
for Weyl spinors harmonics on S 3 is 

dim(n,4) = n(n + l). (2.28) 

In the case of a spacetime with even dimensionality d, the dimensionality of 
the corresponding Clifford algebra of Dirac 7 - matrices is equal to 2 d ; corre- 
spondingly the number of components of Dirac spinors is 2 d / 2 , while the number 
of components of Weyl spinors equals 2( d ~ 2 " 2 . Thus, indices A and A in the 
expression for constant tensor T take 2^ d ~ 2 ^ 2 different values. In such a way the 
problem of calculation of dimensionality of irreducible representation for spinor 
harmonics reduces to the combinatorial counting of the number of possible de- 
composition of the set of n elements into the union of 2 ( - d ~ 2 ^ 2 subsets. One has 

r ( n + 2( d - 2 )/ 2 ) r (n + 2^' 2 - l) 

dim(n ' d) = r (2 (^)/ 2 )i> + i) x r (2 (^)/ 2 )i» • (2 ' 29) 

3. Quantisation of fermion fields in Hamiltonian 
formalism and cosmological wave function 

Following the Dirac procedure for systems with constraints, one may obtain 
from the action Q2.17D the Hamiltonian 

H n (x, x, y, y) = v{xx + yy) + m(yx + xy), (3.1) 

where x, y, x, y obey the Dirac-brackets relations 

{x,x}* = -i, {y,y}* = -i. (3.2) 
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The total Hamiltonian is the sum of the background Hamiltonian, depending on 
the inflaton scalar field and the fermionic Hamiltonian, representing the sum of 
Hamiltonians ( |3.1|) . The total Hamiltonian vanishes, giving in such a way the 
Hamiltonian constraint: 

H + H f = 0. (3.3) 

In the quantum theory, this Hamiltonian constraint becomes the Wheeler-DeWitt 
equation 

[H Q + H f ]\*) = 0, (3.4) 

where is the wave function of the Universe, and the fermionic Hamiltonian is 
obtained from its classical counterpart by replacing the variables with operators 
satisfying the anticommutation relations 

[x,x]+ = l, [y,y]+ = l. (3.5) 

In holomorphic representation they read PR EO] 



- d-d 

x->—,y^—. (3.6) 

ox oy 

With this choice, the wave function is a function of the background variables and 
the unbarred variables. Then there is the usual operator-ordering ambiguity in 
going from the classical constraint Q3.3|) to the Wheeler-DeWitt equation 



Following p2| , |18J , we shall adopt the Weyl ordering, which in this case involves 
the substitution 

With the above choice of operator ordering, the Hamiltonian (|3.1| ) becomes 
the operator 

f d d\ ( d 2 \ 
H n = -v + v\ x— + y— )+m[yx+ . (3.8) 

\ Ox Oy J \ oxoy J 

It is well known J3U], |HT|, |T3] that in the semiclassical approximation the wave 
function of the Universe can be represented as a product of the wave function of 
macroscopical background variables (scale factor, inflaton, etc.)[] and an infinite 
number of partial wave functions describing the perturbations of microscopic - 
in our case fermionic - fields, 

V(t, a, <p\x, y) = * (t, a, <p) U Mt\x, y) ■ (3.9) 



-'^We work in the framework of the reduced phase space quantisation in which the wave 
function explicitly depends on time. 
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These partial wave functions satisfy Schrodinger equations with respect to the 
semiclassical time parameter, 

iip n (t\x,y) = H n if) n (t\x,y). (3.10) 



Before solving ( ft.lC ) with the Hamiltonian (|3.8|) , one should introduce the inner 



product for fermionic wave functions appropriate to the holomorphic representa- 



tion, defined for any pair of functions /, g of the Grassmann variables x, y p8l |29 



(/, 9) = / /0, y)9(x, y)e- x ~ x - y Uxdxdydy. (3.11) 



Integration over x, y, x, y is performed according to usual rules of Berezin inte- 
gration 



dx = 0, / xdx = 1, / dx = 0, / xdx = 1, (3-12) 

and likewise for y and y. 

We are now in a position to look for the solution of the Schrodinger equations 
( |3.10| ). Taking into account the Grassmannian nature of the variables x and y, 



one can easily guess that this wave function is a second-order polynomial in these 
variables. Knowing that the no-boundary wave function of the Universe |Ij is 
given by a Gaussian exponential, we shall look for the solution of ( |3.10|) in the 
following form: 

ip n (t\x,y) = v n + v n xy. (3.13) 



Substituting (|3T3|) and (gj) into (|3~ToD , one can get the following system of 



relations for the functions v n and v n : 



iv n + vv n + mv n = 0, (3-14) 
iv n — vv n + mv n = 0. (3.15) 



These equations obviously coincide with the Dirac equations ( 2.21|) and ( P-20| ) for 



y and x. Now the wave function ( 3.13Q can be rewritten in terms of the first of 



these two functions, the second being expressed in terms of the first one: 

ip n (t\x,y) = v n - — + — xy . (3.16) 

m 

Eliminating v n allows one to write a second-order equation for the basis function 

V n , 

i! n + (-iu + m 2 + v 2 )v n = 0, (3.17) 

which obviously has two independent solutions. The choice of a particular solution 
can be done on the basis of the path-integral derivation of the no-boundary wave 
function which automatically picks out a particular vacuum state for quantum 
fields on the DeSitter background. This will be done in the next section. 
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4. Path integral, Euclidean effective action and 
density matrix of fermions 

The definition of the no-boundary wave function of fermionic fields in the 
holomorphic representation of the previous section can be given by the path inte- 
gral over the fermionic fields in the spacetime ball with the Euclidean spacetime 
metric 

ds% = dr 2 + a 2 (r)dQ 2 3 , a(r) = sin Ht, (4.1) 

H 

which can be obtained from the Lorentzian metric by the analytic continuation 
in the complex plane of time, r = ir/2H + it. For the partial wave functions this 
integral reads 

^(r\x,y) = J Dx Dx Dy Dy e- In[x ^ y >y\ (4.2) 

where I n [x, x, y, y] is the Euclidean action following from fl2.17| ) by Wick rotation 
and by omitting the initial surface boundary term which vanishes because in 
the Hartle-Hawking prescription this surface shrinks to a point satisfying the 
requirement of regularity of all fields. The integration runs over fields with fixed 
values of the unbarred fermionic variables (x, y) at the final (outer) boundary Sf 
labelled by the Euclidean time r. The Gaussian integration gives the answer 

d 2 t 
^(r\x,y) = Det ( - — + v' + v 2 + m 2 ) e~^\ 

= \ W)x[r) + y[r)y[r)\ , (4.3) 

where the prime denotes derivation with respect to Euclidean time, and J^ xt is 
the value of the Euclidean action at the classical extremal, and the functional 
determinant represents the one-loop contribution of quantum fluctuations near 
this extremal. 

The classical Euclidean equations of motion easily follow from their Lorentzian 
counterpart 

—x' + vx — my = 0, (4-4) 
—x! — vx + my = 0, (4-5) 
—y' + vy + mx = 0, (4.6) 
—y' — vy — mx = 0. (4-7) 

The solution of this system of equations in the Euclidean ball < f < r with 
fixed values x{r) = x, y(r) = y at the boundary, 

u(t) u(t) 
m = y ' ( / > y, m = — 1 ' / J x, (4.8) 
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can be given in terms of one function u(t) satisfying the Euclidean version of 

u" + [y' + m 2 + v 2 )u = 0, u(0) = reg, (4.9) 

and the condition of regularity at r = 0. Thus, with this solution the exponential 
of the Euclidean wave function reads — J° xt = xy{u'{r) — uu{r)) / mu{r) . The 
preexponential functional determinant can also be expressed in terms of this 



Euclidean basis function according to the technique of ||, |34|, |35], [15|| . As shown 



in j33], [T^| this determinant, when written as a product of eigenvalues A of the 



differential operator, can be identically transformed as 
d 2 

Det(- — + z/ + z/ 2 + m 2 ) =JJA = E(0), (4.10) 

A 

where E(0) is a particular value at A = of the left-hand side of the equation 

E(X) = (4.11) 

for the eigenvalues of the operator in question. In our case we have a Dirichlet 
problem (zero-boundary condition at the boundary of the Euclidean ball r) for the 
second order differential operator with the A = eigenfunction u(t). Therefore, 
E(0) = u{t) in the equations above, so that the needed functional determinant 
coincides with the basis function taken at the boundary of the Euclidean ball, 

d 2 

Det(- — + v' + v 2 + m 2 ) =u(t). (4.12) 



dr 2 

The resulting Euclidean no-boundary wavefunction of fermions reads 

iv/ i \ (u' — vu\ . . .xy /Atr,\ 

ip {t\x, y) = u exp xy = u + [u — vu) — , (4-13) 

\ mu J m 

and, when analytically continued to the complex plane of time r = n/2H + it, 
gives exactly the wave function (|3.16|) that was obtained above by directly solving 
the Schrodinger equation. The corresponding Lorentzian basis functions 

v(t) = u(n/2H + it), v*(t) = u(n/2H-it) (4.14) 

turn out to be, respectively, the negative and positive frequency basis function of 



the DeSitter invariant vacuum of fermions |36|. This concludes the justification of 
the Euclidean DeSitter invariant vacuum from the no-boundary path integral for 
fermions. Such a justification in the tree- level approximation (disregarding the 



preexponential factors) was achieved for bosons in [31| and for fermions in p2 
Here it is now proven up to the one-loop approximation. Note that for the case 
of fermions, this one-loop contribution brings a nontrivial overall time- dependent 



factor into the wave function ( |4.13| ) which was not correctly considered in f22||p|. 



2 In [p2| the fermionic Hamiltonian was supplied by a complex vacuum energy term, allegedly 
following from the full Whcclcr-DcWitt equation. This certainly leads to nonconservation of 



the norm of the vacuum state differing from ( 3. IS ) by a time-dependent multiplier. Rigorous 
derivation of the one- loop cosmological wave function, made for bosons in [Q, [| [j| , cannot lead 
to such violations of unitarity. 
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Using Grassmann integration rules (|3.12p one can find the positive definite 
square of the norm of the wave function ( |3.16| ): 



Nl = (^„,Vn) = — (\uv n + iv n \ 2 + m 2 \v n \ 2 ) . (4.15) 

Rewritten in terms of barred and unbarred components of the basis function, 
N 2 = \v n (t)\ 2 + \v n (t)\ 2 , v n (t) = u n (n/2H + it), this norm obviously features 
the conservation in time as a consequence of the Dirac equations of motion in 
Lorentzian spacetime. 

Let us now go over to the reduced density matrix of the Universe induced 
by fermions. Integration over fermionic variables as environmental degrees of 
freedom leads for the full cosmological wave function ( |3.9j ) to the expression 

p(a,<p\a\v')=yo%*(UN n j mN^D(a,<p\a',<p'), (4.16) 

where the nontrivial decoherence effects due to fermions are contained in the 
infinite products of norms and the decoherence factor 

D{a,<p\a',<p') = ]lD n {a,<p\a',<p'), (4.17) 



D n (a,tp\a ,(p)- 



m 2 N n N' n 



m 2 + vv' + 





+ zV \-\nv n ) , (4.18) 



where the primes obviously denote that the corresponding quantities are calcu- 
lated on the background of primed macroscopic variables (a', if'). 

The fermionic contribution to diagonal elements of the density matrix reduces 
to the products of the norms of the partial wave functions. Let us show that, 



as in the bosonic case ML E, E, |15|, this contribution is given by the Euclidean 



effective action of fermion fields calculated on the spherical DeSitter instanton: 

1 
2 



niV n = exp(4r), (4.19) 



exp( 



-r) = J] f DxDxDyDye~ In[x >*> y 'y ] . (4.20) 



Here I n [x, x, y, y] is again the Euclidean action, but now it is taken on the whole 
sphere of the DeSitter instanton < r < ir/H on the class of fields regular at 
both poles r = and r = tt/H. Due to the properties of spatial harmonics, this 
regularity implies that the variables simply vanish there, (x, y){r±) = 0, r± = 
0,n/H] that is why both of the surface terms corresponding to these poles are 
also vanishing. The partial contribution to the effective action is obviously given 
by the functional determinant of the same operator as above, 

d 2 

exp(-r n ) = Det(- — + is' + v 2 + m 2 )=l[\ = E n (0), (4.21) 
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but now with the other boundary conditions reflecting the regularity of fields at 
t±. Therefore, the left-hand side of the equation for the zero eigenvalue £7(0) is 
here given by a qualitatively different expression. Let us derive it and, in this 
way, demonstrate that it is given by £7(0) = N 2 - the square of the norm of the 
partial vacuum wave function with the quantum number n, see ( f4.15|) . 

For this purpose, introduce the full Euclidean DeSitter sphere by extending 
the metric ( |4.1|) up to the second pole r = ir/H and smoothly continue the 
operator F = —d 2 /dr 2 + (m 2 + v 2 + v') to the second hemisphere by observing 
that p{t) = v(tt/H — r). Then introduce two sets of basis functions of this 
operator (solutions of the corresponding homogeneous equation Q4.9Q ) 



u'± — (v' + m 2 + v 2 )u± = 0, 
w+(0) = reg, U-(n/H) = reg. 



(4.22) 
(4.23) 



One certainly identifies u + (t) with the Euclidean basis function u(t) introduced 
above and also tries as a candidate for U-(t) the function U\{r) = u(tt/H — r) 
used in the bosonic case in flT5| . Now, however, this turns out to be impossible, 
because in view of the relation i/(tt/H—t) = —v'(t) this function satisfies another 
equation, u'[ — (— v' + m 2 + v 2 )u\ = 0. The remedy is to take first the new function 
U2 = (u' — vu) j m 2 which, as can easily be verified, satisfies this equation with 
the negative sign of u', v!^ — (— v' + m 2 + v 2 )u 2 = 0, and then construct its mirror 
image with respect to the equatorial section of r = n/2H, U-(r) = 112(11: /H — r). 
Then this function will satisfy the original equation and will be, by construction, 
regular at the pole r = ir/H. Thus the set of two basis functions of our operator 
regular at the opposite poles of the DeSitter instanton is given by 



u+ (t) = u(t), u_(r) 



u 



vu 



77T 



(4.24) 



In terms of these two basis functions, the equation of the zero eigenvalue of the 
operator can be written as the linear dependence of two two-dimensional vectors 
composed of these functions and their derivatives, 



u + 


(t)~ 






(r)~ 


u> + 


(r)_ 


— o 


u'_ 


(r). 



0. 



The left-hand side of this equation can also be written as 
E(0) = (u-u' + — u + u'_^(t). 



(4.25) 



(4.26) 



Note that in view of the Wronskian relation for the operator in question, this 
quantity is independent of r. Therefore, when calculated at the equator r = 
ir/2H, on using the equation for u(r) and its corollary u'_{r) = (l/m 2 )(—mu + 
(u' — vu)v){tt / H — r), the quantity E(0) reads 



£7(0) = — r \m 2 u 2 + (u' 



vu] 



t=tt/2H 



(4.27) 
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In view of the analytic continuation to the Lorentzian time r = ir/2H + it, (4.14) 



at r = tt/2H - the point of "nucleation" of the Lorentzian spacetime from the 
Euclidean hemisphere - we have u(n/2H) = v(0) = v*(0), u'(n/2H) = —iv(0) = 
iv*(0) and 

E(0) = N 2 (4.28) 



with the square of the norm given by the expression (|4.15 ) - the quantity which 



is in its turn time-independent in view of the Dirac equation for basis functions 
in Lorentzian spacetime. 

This finally proves the relations (|4.19|) -( |4.20 ) and brings the fermionic contri- 



bution to the reduced density matrix in quantum cosmology to the form 

PfermionK fWf') = eX P ~ D ^ ^K' <p')i ( 4 - 29 ) 

analogous to the bosonic contribution obtained in our previous paper 

This concludes the proof of the Euclidean effective-action algorithm for the 
diagonal element of the density matrix - the distribution function of inflationary 
cosmological models - derived in the bosonic case in the series of publications 
@, H [| , and now extended to fermions. For the off-diagonal elements the effective 
action factors exp(— ^r) and exp(— \r') also modify the usual expression |[S 



but these factors are constant in t and, therefore, become negligible for late 
Lorentzian times, provided the decoherence effects due to D(a, <p\a', <p') rapidly 
grow with time. Therefore, in what follows we shall concentrate on the time 
evolution of the decoherence factor D(a,ip\a',<p'). 



5. Decoherence due to massless versus massive 
fermions 

The structure of the decoherence factor for fermions is much more complicated 
than the corresponding one for bosons (cf. |15|). Their quantum state ( |3.16| ) 
generated according to the no-boundary prescription can be interpreted as the 
vacuum state corresponding to the following creation and annihilation operators: 

\uv n - iv n \ d \v n \m 



\J\w n 




2 _|_ m 2 


V n \ 2 




\w n - 


~ iv n \ 






-iv n 


2 _|_ m 2 


\ V n\ 2 




vv n - 


- iv n \ 








2 _|_ m 2 


V \ 2( 

u n \ 




\vv n - 


" iv n \ 





^\vv n - iv n \ 2 + m 2 \v n \ 2 
\v n \m d 
vv n — iv n \ 2 + fn 2 \v n \ 2 ® x 



y 



,..„.. „ d \v n \m 
b n = i. = = = — ?. V> 

I 2 I /yyi 2 Li I 2 



y/\uv n - iv n \ 2 + m 2 \v n \ 
sj\vv n - iv n \ 2 + m 2 \v n \ 2 \j\vv n - iv n \ 2 + m 2 \v n \ 2 °V 



14 



In the massless limit m = 0, these operators have a very simple form: 
d_ 

6„ = (5.2) 

and the vacuum state contains only the Grassmann-independent component fol- 
lowing from the Schrodinger equation (|3.10|) : 

ip n = exp (i / u(t)dt). (5.3) 



From the structure of the wave function ( |5.3|) one can see that the decoherence 
factor does not have a real part and hence, decoherence effects are absent. This 
phenomenon was already described in It is not surprising because the ab- 
sence of decoherence was observed not only for conformally invariant massless 
fermions but for all conformally invariant fields such as the electromagnetic field 



and conformally-coupled massless scalar field |L5 . 

Let us turn now to the case of a massive spinor field. To get the Hartle- 
Hawking vacuum one should take the wave function (|3.16| ), where the function 



v n can be obtained by analytic continuation from the solution of the Euclidean 
counterpart of the second-order equation ( |3.17| ). One should choose a solution of 



the Euclidean equation which is regular at the Euclidean hemisphere and serves 
as a germ for the Lorentzian Universe |5|, ||. Such a solution can be found for 
example in |36| and has the following form: 

v n = (1 + % sinh Ht) n/2+3/ \l - i sinh Ht)- n/2+1/4 

x 2 Fi(l + im/H, 1 - im/H; n + 2; (1 + i sinh#t)/2). (5.4) 

We are interested in the case of large masses because at the beginning of inflation 
these large masses can be generated by Yukawa interactions between fermions 
and the inflaton scalar field (see, for example ||, 0). In this case it would be 
convenient to have a uniform asymptotic expansion for the solution of ( |3. 17] ). 



Unfortunately, in contrast to Legendre or Bessel equations (cf. |J7|, [33|]), such 
an expansion cannot be written down in closed form. However, one can write 
down the main terms of this expansions responsible for the leading structure in 
the decoherence factor ( |4.18| ). The leading- in- mass m part of the basis-function 
frequency looks like 
d , i 



— In v n ~ —V n 2 + m 2 a 2 . (5-5) 
at a 

Substituting it into the formula for the decoherence factor D n , one gets the fol- 
lowing expression for the leading term of the amplitude for the total decoherence 
factor: 

\D(a, <p\af, v0| = II \Dn\ = exp (- ^"^ /) , (5.6) 
- n(n + l) 

^Srw (5 - 7) 
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It is interesting to notice that in the opposite limit of small masses m one can 
immediately use the hypergeometrical expansion ( |5.4j ), and the leading term in 
mass again has the structure (|5.6|) . This expression for small mass was obtained 
earlier in 



It is clear that the expression ( |5.6|) contains linear ultraviolet divergences. 
It is in sharp contrast with the case of boson fields, studied in our preceding 
paper [|nj. Indeed, in the case of boson fields, using conformal parametrisation 
of fields, it was possible to avoid ultraviolet divergences (see also |T6j] ). Now, in 
spite of using the conformal parametrisation, this is not the case here. Moreover, 
we cannot freely redefine the fermion field by multiplying it with a power of the 
cosmological factor a, because the measure on the corresponding Hilbert space of 
states (see (|3.11| )) then starts depending on this multiplier. With two scale factors 
a and a' related to the arguments of the density matrix this would mean having 
two different integration measures, neither of which can be used for constructing 
the reduced density matrix. Thus, we should try to renormalise the decoherence 
factor by some regularisation procedure. This will be done in the next section by 
means of dimensional regularisation [P 



6. Renormalised fermionic decoherence factor 

We would like to calculate now the renormalised expression for the decoher- 



ence factor (|5.6j ). This calculation is reduced to the calculation of the regularised 
sum ( |5.7| ). If this sum is positive, then there is a decoherence effect, while a 
negative sign of this quantity would imply intrinsic inconsistency of the obtained 
density matrix which is likely to become an unbounded operator ||15|| . 

In the framework of dimensional regularisation, dim(n, 4) = n(n+l) should be 
substituted by dim(n, d) from ( |2.29|) , and the eigenvalue of the three-dimensional 
projection of the Dirac operator, {n + 1/2), should be substituted by that of the 
(d — 1) - dimensional Dirac operator, (n + (d— 3)/2). Thus, one has 

^ dim(n,4) 

1 = £ (»+(«■- d/2)»- {6A) 

To calculate the sum of this type, one can use the Euler-Maclaurin formula: 

I = h + h, (6.2) 
r°° dimfa, d) 

h = l dy (y + (d-Wr (6 - 3) 

. f°° dy ( dimfa, d) dim(-iy, d) \ , , 

2 l Jo exp(27q/)-l \(iy + (d - Z) / 2) 2 {-iy + (d - 3)/2) 2 J " K ' } 

The last integral (see (|6.4|) ) is convergent at d = 4 and equals pf] : 

1 r* dyy n * - 1. (e.s) 



2 Jo (exp(27ry) - l)(y 2 + 1/4) 2 8 
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The integral I\ in ( |6.3j ) is linearly divergent and for its regularisation it is 
necessary to use the analytic-continuation procedure and integration by parts 
typical for dimensional regularisation |20|, [15| . Now it is convenient to make a 
change of variables from y to x = 1/y and to introduce the new quantity 

_ 4dim(l/x, d) 2d/2 _ 2 
A^j [x(rf _ 3) + 2] 2 X > ^-bj 

which is finite at x — > oo. Now, the integral I\ can be rewritten as follows: 

Ji = dxf(x,d)x~ i2d/2 - 2 \ (6.7) 

Integrating two times by part and discarding the terms corresponding to power- 
law divergences, one gets 

h = ( 2^-3) 1 (2^-4 ) r^ , ^-''^- |2J/2 - 4 ' < 6 - 8 > 



Then, expanding expression (|6.9| ) in powers of (<i — 4) one has the following 
regularised expression: 

h = ~[ dx In xf"(x, 4) - + /'(0, 4) - 4M, (6.9) 

where , d denotes the derivative with respect to the dimensionality of spacetime. 
The direct calculation shows that 

f(0,4) = 0, (6.10) 

which corresponds to the absence of logarithmic ultraviolet divergences. Then 

1 

dx\nxf"(x,A) = +-, (6.11) 

and the mixed derivative is given by 

/; rf (0,4) = -l + 21n2. (6.12) 
Substituting (|6U01) - (|dU21) into (gj), one has 
l-31n2 

h = , (6.13) 

and combining ( |6.13| ) with ( |6.5| ) one finally has 
l-51n2 7T 2 

/ = + — < 0. 6.14 

21n2 8 v ' 

Thus, we have found a negative result which would be in strong contradiction with 
basic properties of a density matrix, since Trp 2 would diverge. The next section 
will show how the divergences can be remedied without leading to inconsistencies. 
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7. Bogoliubov transformation for fermions and 
decoherence 



We have already seen in Sec. 5 that ultraviolet divergences arise in the 
decoherence factor of the reduced density matrix coming from fermionic pertur- 
bations, and in Sec. 6 we have obtained a dimensionally-renormalised expression 
which manifested wrong, anti-decoherent behaviour, breaking such fundamental 
properties of the density matrix as the normalisabity of its square which is equiv- 
alent to its positive-defmiteness and boundedness. It happens in spite of working 
from the beginning in the conformal parametrisation for the spinor field, which 
was efficient in eliminating ultraviolet divergences in decoherence factors induced 



by bosonic perturbations [T§, [15[ . Moreover, we do not have the possibility to 
change the parametrisation of fermion fields by simply multiplying them by some 
power of the scale factor a, as it was done for bosons [|15j . This happens because, 



as was mentioned in Sec. 5, the inner product of wave functions depending on 
fermion fields depends on these fields. 

However, we still have some additional freedom of redefining the environmen- 
tal variables by performing some kind of Bogoliubov transformations. This will 
lead to well-defined expressions. At the end of this section it will be shown that 
this transformation is analogous to the well-known Foldy-Wouthuysen transfor- 



mation 21 



Thus, let us make the following Bogoliubov transformation 

x = ax + (3y, 

y = iy + Sx, 

x = ot*x + P*y, 

y = j*y + 6*x. (7.1) 



Substituting the expressions (|7. 1|) into the action ( [2.17D and requiring the conser- 



vation of the form of the kinetic term xx + yy, one has the following requirements 
for the Bogoliubov coefficients a, j3, 7, 6 : 

a */3 + 7 *5 = 0, 
\a\ 2 + \f3\ 2 = 1, 

| T | 2 + |<5| 2 =1. (7.2) 
It is convenient to choose the Bogoliubov coefficients in the form: 

a = 7, 

P = -S. (7.3) 



Now, substituting the transformed Grassmann variables ( |7. 1|) into the other 
terms of the action ( |2.17| ), taking into account ( [7. 2| ) and ( |7.3| ) and arriving at the 
Hamiltonian formalism, one can get the new Hamiltonian 

H n (x, x, y, y) = v{xx + yy) + frfyx + mxy, (7.4) 
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where 



m 



i{aa* + (3(3*) + v(\a\ 2 - \(3\ 2 ) - m{a(3 + a*[3*), v* = v 
-- m(a* 2 - (3 2 ) - 2va*(3 + i((3a* - pa*). 



(7.5) 
(7.6) 



Thus, we have found a new Hamiltonian which has the same form as the old 
one (|3.1|) , with the only difference that v and m are replaced by new functions v 
and rh depending on Bogoliubov coefficients. Besides, the "mass" rh is complex 
now (p remains real as one can check by using the relation ( |7.2| ) for Bogoliubov 
coefficients). 

The Dirac equations now read 



ix + vx — my = 0, 
ix — vx + m*y = 0, 
iy + vy + rhx = 0, 



iy — vy — m*x = 0. 



(7.7) 
(7.8) 
(7.9) 
(7.10) 



For generic Bogoliubov coefficients (a n , (3 n ) depending on the quantum num- 
ber n, the decoherence factor (|4.18j) becomes 



D(a, p\a f Lp') = Yl D n (a, ip\a', tp'), 
D n (a, ip\a'(p' 



(7.11) 



VnV'* 



m n m' n N n N' n 



rhjht 



vv + 



' d_ 
dt 



In v. 




—iv 



In v'* + iv' 



' d_ 
dt 



lnv n 



(7.12) 



with the new mass parameter m n also depending on n (v was n-dependent even 
in the old parametrisation), while the approximate expression for a potentially 
divergent part of this factor (|5.6|) becomes 



\D(a,ip\a',tp')\ = n \D n \ = exp 



[a — a 



A2 



E 

n=l 



n(n + l)m n m' r 



in 



l/2) s 



(7.13) 



replacing ( p.6[ ) 

The dependence of the new functions v and m on Bogoliubov coefficients gives 
new opportunities regarding the finiteness and consistency of the decoherence 
factor. In particular, the modified mass rh can be demanded to vanish. This 
would, however, lead to an absence of decoherence. For this reason, a physical 
principles must be invoked to select an appropriate Bogoliubov transformation: 
We demand that there be no effect for a stationary spacetime, in accordance 
with the idea [15| that decoherence is connected with particle creation. The 
decoherence factor should thus contain terms that vanish when a vanishes. Since 
one the main features of decoherence is its irreversibility § - and the escape of 
created particles is, similar to usual scattering situations, an irreversible process 
- we think that the above principle captures the physical effect of decoherence. 
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We therefore choose real a and ft satisfying the equation (this is the (at,ft)- 
independent part of in) 

m(a 2 n -ftl)-2ua n ft n = 0. (7.14) 

One finds the following solutions: 

a n = ( ^+^±A 1/2 f (7 . 15 ) 

Pn={^S^-T- (7-16) 



2^ 



m A + v 



If v is independent of time (a stationary Universe), this choice of Bogoliubov 
coefficients results in a vanishing m; however, for a nonstationary Universe there 
is an imaginary contribution proportional to time derivatives of Bogoliubov co- 
efficients, 

. _ m 2 i> ( Vm 2 + v 2 + v \ 1/2 (7 -\7\ 

an ~ 4(m 2 + 1/2)3/2 I 2Vm 2 + v 2 ) ' 1 ' 



m 2 v 



-1/2 



o _ — ( Vrn 2 + u 2 -u " 

Pn ~ A{m 2 + v 2 fl 2 \ 2v^T^ J ' 1 1 

Therefore, 

mi> 

™ n = i ( n— — 2T, (7.19) 
2(m^ + v l ) 

and at large values of n this mass parameter behaves as 

m n ~ - . (7.20) 

n 

Then, turning to the expression (|7.13| ) one can see that the dependence of 
fh on n presented in ( [7.1 9| ) provides the suppression of ultraviolet divergences in 
the modified expression (|7.13|) . The terms in the sum in the exponent of ( [7.1 3|) 
approach for large n the expression 

n(n + l)m n m* n m 2 a 2 

{n + l/2) 2 ~* 4n 2 ' 1 ' j 

and the sum is explicitly convergent. 

We want to conclude this section by comparing the above Bogoliubov transfor- 
mations with their counterparts in a stationary flat spacetime. As it was already 
mentioned for the case of stationary background, coefficients a and ft from Eqs. 
( 7.16[ )-( [7.16 ) eliminate mass exactly. The form of these coefficients reminds that 
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of a Foldy-Wouthuysen transformation PTJ, but really their role is just the oppo 



site. To explain it, let us remember that in the canonical basis (cf. with (|2.9|) in 
Sec. 2), a Dirac spinor is represented in the form 

*c = g) , (7.22) 
where two-component spinors satisfy the following Dirac equations: 

~ M A 

% ~dt = apx + m( P' 

To disentangle these equations and to get rid of the so-called odd terms propor- 
tional to ap, one can use a Foldy-Wouthuysen transformation 

s 1 Fn ., (7.24) 



\P\ 



a: 





\P\ 



r-^z- (7-25) 

This transformation also digonalises the Hamiltonian which is now represented 
by a direct sum of two non-local Hamiltonians ±y/m 2 + p 2 : 

H = f3^/m 2 + p 2 . (7.26) 

It is clear that these square roots cannot be represented in configuration space 
by a finite set of differential operators. It is obvious also that our modified 
Hamiltonians ( [7.4[ ) are also combined into a non-local Hamiltonian because all 
parameters u, rh of these partial Hamiltonians depend on the wave number n. 
However, as a matter of fact, we are working in the Weyl representation for 



Dirac spinors ( |2.8| ) and 7 matrices ( |2.7|) . In this representation, especially, massive 
terms are mixing in the Dirac equations ( p,12j ) and our task is to get rid of them. 
This is exactly the problem that was studied in this section in two-component 
formalism; it can thus be called "anti- Foldy-Wouthuysen" transformation and can 
be in a traditional (four-spinor) formalism presented by the following operators: 

* = e T ^, (7.27) 



where 



T= ~ Q og 




■x 1 Tn 

9 = — arcsin = =. 7.28) 
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Correspondingly, the Hamiltonian is turned by this transformation into 

m 

and is again apparently non-local. 



8. Conclusions 

Let us briefly recapitulate the main results of the present paper. We devel- 
oped the formalism of the one-loop no-boundary wave function of the Universe 
in a cosmological model with fermions establishing the connection between the 
path integral, effective action and semiclasssical Schrodinger equation. It is worth 
adding here that these results are valid for the tunneling wave function || as well 
due to the fact that the structure of the microscopic perturbation part of these 



two wave functions is the same |39], |40j . We studied the expression for the reduced 
density matrix and showed that the amplitude of its decoherence factor reveals 
ultraviolet divergences, which was previously observed in The renormalised 
expression for this decoherence factor which was obtained by dimensional reg- 
ularization has a wrong sign and violates the properties of density matrix as a 
bounded operator. However, by making special Bogoliubov transformation for 
fermionic variables one can get a finite expression for this decoherence factor. 



Thus, as in the case of bosonic variables |15| , the definition of the degrees of 
freedom constituting an environment is of crucial importance for the construc- 
tion of consistent expressions for the reduced density matrix (see also |]16 , [17|). 



However, there is an essential difference between the fermionic and bosonic cases. 
In the case of bosonic degrees of freedom we have a well-defined criterion for the 
choice of the demarcation line between the system under consideration and the 
environment. Unfortunately, the situation is not so clear for fermionic degrees of 
freedom. We have put forward the suggestion to fix the Bogoliubov transforma- 
tion by demanding that the corresponding decoherence effect would be absent if 
the spacetime were stationary. 

One problem is that any choice of this transformation breaks locality of a 
fermion field; this can create difficulties at the consideration of fermions intere- 
acting with other fields. However, in contrast with bosons when the logarithm of 
a decoherence factor is proportional to m 3 in the limit of large mass, for fermions 
that is proportional only to m 2 . This fact gives support to the conclusion that 
bosons are much more efficient than fermions for the quantum-to-classical tran- 



sition 18 



In summary, we have extended the Euclidean effective action algorithm for the 
quantum distribution function of cosmological models and their density matrix 
to fermions; this was previously only known for bosons j|, [5], [], [RJ. Together 



with |15[ this gives a complete account of the quantum-to-classical transition at 



the onset of inflation - the beginning of our classical Universe. 
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